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ABSTRACT 

Gamma Ray Bursts are detectable in the y-ray band if their jets are oriented towards the observer. However, for each GRB with a 
typical j), there should be ~ 2/# 2 ct bursts whose emission cone is oriented elsewhere in space. These off-axis bursts can be eventually 
detected when, due to the deceleration of their relativistic jets, the beaming angle becomes comparable to the viewing angle. Orphan 
Afterglows (OA) should outnumber the current population of bursts detected in the y-ray band even if they have not been conclusively 
observed so far at any frequency. We compute the expected flux of the population of orphan afterglows in the mm, optical and X-ray 
bands through a population synthesis code of GRBs and the standard afterglow emission model. We estimate the detection rate of OA 
by on-going and forthcoming surveys. The average duration of OA as transients above a given limiting flux is derived and described 
with analytical expressions: in general OA should appear as daily transients in optical surveys and as monthly/yearly transients in the 
mm/radio band. We find that ~ 2 OA yr -1 could already be detected by Gaia and up to 20 OA yr -1 could be observed by the ZTF 
survey. A larger number of 50 OA yr~* should be detected by LSST in the optical band. For the X-ray band, ~ 26 OA yr -1 could be 
detected by the eROSITA. For the large population of OA detectable by LSST, the X-ray and optical follow up of the light curve (for 
the brightest cases) and/or the extensive follow up of their emission in the mm and radio band could be the key to disentangle their 
GRB nature from other extragalactic transients of comparable flux density. 

Key words, stars: gamma-ray bursts: general, relativistic processes 


1. Introduction 


Gamma Ray Bursts (GRBs) are cosmological sources which 
signpost the birth of stellar mass black holes or the most intense 
magnetic fields harbored by compact objects (magnetars). They 
are detected as short, highly variable, flashes of y-rays (prompt 
emission) with duration that is typically in the range 0.1 to 1000 
seconds followed by a smoothly decaying long lived emission 
at X-ray, optical, radio frequencies (the afterglow - Costa et al. 
19971: van Paradijs et al. ll997l) . In the standard fireball model the 
prompt emission is interpreted as due to internal dissipation (ei¬ 
ther through relativistic shocks or magnetic reconnection) while 
the afterglow is produced by the deceleration of the relativistic 
outflow by the circum burst interstellar medium. 

Theoretical arguments (e.g^ Sari & Pi ra n 1999i) and direct 
observational evidences (e.g. iMolinari et alT12007 ) suggest that 
the outflow of GRBs is relativistic with typical bulk Lorentz fac¬ 
tors To~ 10 2 3 . The most luminous and energetic GRBs seem to 
have larger To (e.g. iLiang et aDl2010t ?). Another key property 
of GRBs is the presence of a jet. Invoked to reduce the other¬ 
wise huge isotropic equivalent energies by a factor proportional 
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to 6 2 (i.e. the jet half opening angle), jet angles have been esti- 

jet 

mated in a few tens of GRBs from the steepening of the optical 
and/or X-ray light curve a few days after the prompt emission. 
This steepening is interpreted as the time when, due to the de¬ 
celeration of the outflow, the relativistic beaming oc 1/T equals 
the geometric collimation f> . Modelling of the outflow dynamics 
(iBlandford & McKeeIl976l) allows us to infer the GRB jet open- 
ing angle b dRhoadsI 1 1 999l> , The collimation c orrected energies 
were found to be distributed around 10 51 erg IlFrail et alil200lh 
with a smaller dispersion with respect to that of the isotropic 
equivalent energies (but see lGehrels et al.[|2009l) . 

The satellites/instruments deputed to the detection of GRBs 
observed hundreds/thousands of GRBs at an average rate of ~ 
0.3 day -1 . However, if GRBs have a jet with the energy and bulk 
Lorentz factors which are constant within the jet as a function 
of the angle from its axis and Tq (pointed radially as the expan¬ 
sion of the outflow within the jet) is relatively large (typically 
a few hundreds), we can detect only those bursts whose jet is 
pointing at the Earth. Indeed, since the highly relativistic mo¬ 
tion results in a strong forward beaming of the emitted radia¬ 
tion, the flux directed at the Earth is dramatically reduced when 
#view>£> (where 0 v j ew is the viewing angle between the jet axis 
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and the line of sight). These events (which are the most numer¬ 
ous due to the jet orientation probability being oc sin 0 v i ew ) go 
undetected as prompt y-ray bursts. However, during the after¬ 
glow the bulk Lorentz factor decreases with time (as the outflow 
is decelerated by the external medium). There will be a charac¬ 
teristic timescale when the relativistic beaming oc 1 /T equals the 
observer viewing angle 6 V j ew and the (afterglow) radiation can be 
seen. These events, missing the prompt emission but detected as 
afterglows, are called orphan afterglows (OA - hereafter) and, 
for typical opening angles of GRBs of a few degrees, e.g. [t~0.1 
rad, they should outnumber the population of GRBs (by a factor 
OC (1 - COS0jet) _1 ~ 200). 

Therefore, OA should be detected as transients but their 
association with GRBs is made difficult due to the lack 
of a prompt high-energy emission. Despite specific stud- 


ies have been searching for OA in X-r 

ay su 

sun 

rveys (iGrindlav 

1999 

; Greiner et al. 2000) in optical 

/eys (vreeswiik 

2002 

; Rau et al. 2006; Malacrino et al. 20071 

Rau et al. 2007| 

and in the radio band (Levinson et al. 

2002 

Gal-Yam et al. 

2006; Bannister et al. 2011; Bell et al. 

2011 

Bower & Saul 

2011; Croft et al. 2010; Frail et al. 2012; Carilli et al. 2003; 

Matsumura et al. 2009; Lazio et al. 2010j), 

no OA has been con- 


clusively detected so far. 

Non detection s of OA are in agreement with current theo¬ 
retical predictio ns jTbtani& ^anaifescu|2002 , Nakar et al. 120021: 
IZou et al.ll2007t iRossi et aTT 20081 iMetzger et al.ll2015l) . How- 
ever, these works either extrapolated the properties of a few 
know n GRB afterglows to the orphans (e.g. Totani & Panaitescu 
l2002t> or assumed basic prescriptions for the known GRB pop¬ 
ulation properties or for the afterglow emission model. We have 
recently developed a population synthesis code dGhirlanda et al.1 
1201 3 al) which, cou pled with the most d etailed model for the af¬ 
terglow emission (Ivan Eerten et al. 112012l). allowed us to predic t 
the properties of the population of OA dGhirlanda et al.ll2014l) 
reproducing a large set of observed properties of the population 
of the “Earth-pointed” GRBs. We have, so far, considered the 
radio band predicting that the Square Kilometer Array (SKA), 
reaching the p}y flux limit, could see up to ~0.2— 1.5 OA deg~ 1 2 
yr -1 dGhirlanda et al .1120 141) . Alternatively, th e non detection of 
OA could be due to the structure of the jet dRossi et al.l [2008: 
ISalafia et al.ll20l5b . 


We are entering the era of large synoptic surveys which will 
monitor large portions of the (if not the whole) sky with unprece¬ 
dented sensitivities. OA are potentially in the list of transients 
that these surveys will detect, but specific predictions on the rate 
depend on the true rate of the population of OA (and also their 
duration) and on the survey characteristics (area of the sky cov¬ 
ered, timescales, limiting flux). Here we derive the flux distribu¬ 
tion of OA in the X-ray, optical and mm band (§3) based on our 
recent population synthesis code (§2). We also study the average 
duration of the population of OA as a function of the survey lim¬ 
iting flux (§3). We compare with current limits of OA in these 
bands and make predictions for on-going and forthcoming sur¬ 
veys (§4). Standard cosmological parameters (h = O a = 0.7) for 
a flat Universe are adopted throughout the paper. 


2. Orphan Afterglow emission 


The only difference between OA and GRBs is the orienta¬ 
tion of their jets with respect to the line of sight. This allows 
us to use all the known properties of GRBs detected so far in 
the y-ray band and with well studied afterglow emission to infer 



L°g(F) [«Jy] 

Fig. 1 . Flux density cumulative distributions in the optical (R - solid 
red line) X-ray (at 3 keV - solid blue line) computed at 11 h after the 
start of the GRB emission for the BAT6 Swift complete sample (adapted 
front Melandri et al. l2014l and D’Avanzo et al. l2012l for the optical and 
X-ra y band, respectiv ely). The radio (at 8.4 GHz - solid green line) is 
from IGhirlanda et afl d2013bl ). also for the BAT6 sample. The results of 
the population synthesis code with p = 2.3, e e = 0.02 and e B = 0.008 
are shown with the dashed lines. The shaded regions represent, for each 
band, the results obtained with (p,e e ,e B ) = (2.3,0.01,0.001) for the 
lower boundary and (p, e e , e B ) = (2.3,0.05,0.01) for the upper bound¬ 
ary. The X-ray and R band fluxes of the Swift BAT6 sample (solid blue 
and red line respectively) have been corrected for absorption (D’Avanzo 
et al. 2012; Campana et al. 2012) and for dust extinction (Melandri et 
al. 2014; Covino et al. 2013), respectively. 


the emission characteristics of OA. In particular, OA are normal 
GRBs with their jet oriented so that d v iew>f>- 

In order to predict the properties of OA we first need a 
model describing the entire population of GRBs distributed 
in the Universe. We use the population code developed re- 
ce ntly in GhiiTanda etaLl (1201 3al G13 hereafter) and extended 
in lGhirlanda et alJ (12014 G14 hereafter), called PSYCHE (Pop¬ 
ulation SYnthesis Code and Hydrodynamic Emission model). 
PSYCHE generates burst s with a redshift z (assigne d following 
the GRB formation rate - iHopkins & Beacomll2008l ). a jet open¬ 
ing angle jj, a bulk Lorentz factor To. The latter two parameters 
are drawn from two log-normal distributions with median val¬ 
ues 5.1° and 90, respectively. These distributions were obtained 
(G13) in order to reproduce: (a) the E p - E IS0 correlation of a 
comp l ete (flux limited) sample of Swift bursts (Salv aterra et al.l 
12012b iNava et aPl2012l) : (b) the flux distribution of GRBs de¬ 
tected by BATSE and GBM/Fermi; (c) the detection rate of 
GRBs by Swift, Fermi and BATSE. The isotropic equivalent en¬ 
ergy £ is0 and rest frame peak energy E p are obtained, once |r and 
Tq are extracted, assuming a universal comoving frame, collima- 
tion corrected, energetic (see G13 for details Q 


Bursts are then assigned a viewing angle 4 v j ew (according to 
the sin0 v i ew probability distribution) representing the orientation 
of the jet with respect to the line of sight. Within the GRB pop¬ 
ulation simulated by PSYCHE, there are bursts that can be de¬ 
tected as y-ray events because 4 v j cw <|t and those that that can be 


1 The distribution of jt and r 0 are derived self-consistently in G13 to 
reproduce the y-ray properties of GRBs. It is shown in G13 that even 

changing the assumed values of the comoving frame energetics by a 
factor of 10 the two inferred distributions are modified accordingly so 
that the population of GRBs has similar characteristics in terms of en¬ 
ergetics and opening angles 
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detected as orphan afterglows because 0 v iew>f>- The latter are the 
subject of this work. 

In general, since for highly relativistic GRBs seen off-axis 
only the afterglow can be seen (§1), we need to simulate for 
each burst its afterglow emission. To thisjiim PSYCHE imple¬ 
ments the numerical code B OX FIT (Ivan Eerten & MacFadvenl 
1201 It Ivan Eerten et alJl2012l VE12 hereafter), which is based 
on numerical 2D simulations of the jet dynamics and assumes 
synchrotron emission from shock accelerated electrons as the ra¬ 
diation mechanisms of the afterglow phase. The VE12 code con¬ 
siders GRBs with a constant circum burst density. In this stan¬ 
dard afterglow model, the external shock emission depends on a 
set of micro-physical parameters: the index p of the energy dis¬ 
tribution of the shock accelerated electrons, the fraction of the 
dissipated energy distributed to electrons e e and to the magnetic 
field 6b . Finally, the value of the density n for the circumburst 
medium, is assigned. These are the free parameters determining 
the afterglow emission of each simulated burst. 

2.1. Setting the micro-physical shock parameters 

Distributions of the micro-physical shock parameters ( p , e e , es) 
are poorly constrained directly from the observations. Dense 
multiwavelength sampling of the afterglow light curve from 
early times to days after the burst explosion is availa ble for a 
limited number of bursts. IPanaitescu & Kumar i ( 2000 ) first de¬ 
rived the values of these parameters through the modelling of 
the afterglow light curves of 10 GRBs in the pre-Swift era. 
The sa me pre-Swift bursts were used to derive the properties 
of OA dTotani & Panaitescull2002h . However, Swift follow up in 
the X-ray band and optical monitoring campaigns have shown 
in the last years that the afterglow emission can be very differ¬ 
ent from burst to burst. The X-ray and optical luminosities at 
0.5 days after the bur st have proven (e.g. iD’Avanzo et alJl2012t 
iMelandri et al.l 120141) to be more dispersed than initially found 
based on a limited number of events. Moreover, a possible differ¬ 
ent origin of the X-ray and optical emission, with the latter being 
more genuinely afterglow, has been proposed ((Ghisellini et all 
120091 ). 

The aim of this work is to estimate the flux level of the pop¬ 
ulation of OA as a whole in order to compare it with the lim¬ 
its of current and future surveys in different bands. To this aim 
we can use average values of the micro-physical parameters (p, 
e e , €b) for the entire simulated GRB population. In order to as¬ 
sign these values, we consider the observed flux of the after¬ 
glow of the GRBs composing the complete Swift sample (BAT6 
- ISalvaterra et alJl2012li . This is a flux limited sample of GRBs 
and it turns out to have a high level of completeness in red- 
shift (97%). Therefore, no bias induced by the redshift estimate 
should be affecting it. Similarly the high flux cut ensures that it 
is free from any threshold bias related to the GRB detector (BAT 
on board Swift in this case). Fig.|T|shows the cumulative distri¬ 
bution of the optical flux (at 12 h after the burst) of the BAT6 
sample (red solid line - adapted from lMelandri et aPl2014l) . the 
cumulative distribution of the X-ray flux (at 3 keV and at_l lh) of 
the BAT6 (blue solid line - adapted from lD’Avanzo et~aPl2012l) 
and the cumulative distribution of the radio flux (at 8.4 GHz be¬ 
tween 1 and 6 days) of the BAT6 (green solid line -adapted from 
iGhirlanda et al.ll2013al ). These are the observed distributions we 
aim to reproduce assigning to the simulated GRB population a 
set of values for the micro-physical parameters (p, e e , 6b). 

Since the BAT6 sample contains GRBs detected by Swift in 
the y-ray band, we work on the simulated population of GRBs 
with $ v iew—|t. We select among them the bright events with the 


Table 1 . Parameters of the linear fits to the average duration of OA 
above flux threshold (Fig. [3j. Fits parameters of the formula: log(< 

T -^days) — q + tn log(F ]im ,mjy) 


Band 

m 

q 

X-ray (3 keV) 

-0.28 

-1.54 

R (7000 A) 

-0.36 

-0.72 

443 GHz 

-0.44 

0.44 


same flux cut adopted for the BAT6, i.e. bursts with a peak flux 
larger than 2.6 ph cm~ 2 s _1 integrated in the 15-150 keV en¬ 
ergy band. We assume that the circumburst density n is dis¬ 
tributed uniformly between 0.1 and 30 cm -3 (also in this case 
the choice is done in order to reproduce the BAT6 optical and 
radio flux distributions - Fig. Q} and assign to each simulated 
burst a value of n randomly extracted from a uniform distribution 
within this range. We keep the other three parameters (p, e e , ef) 
fixed. The fiducial values for these three parameters are obtained 
by reproducing the flux distributions of the BAT6 sample shown 
by the solid lines in Fig. [T] at three different frequencies. With 
p - 2.3, e e = 0.02 and e# = 0.008 (as already discussed in 
Ghirlanda et al. 2013a, 2014) we obtain for the simulated pop¬ 
ulation of GRBs cumulative flux distributions (dashed lines in 
Fig. 0} which nicely match the radio and optical flux distribu¬ 
tions (solid lines in Fig. Q} of the real GRBs of the BAT6 sam¬ 
ple. If we assume smaller/larger values for e e and es we obtain 
smaller/larger fluxes in both bands (solid shaded regions in Fig. 

CD- 

Also shown in Fig. Q] (solid blue line) is the X-ray flux 
c umu lative distribution at 3 keV computed at 11 h (from 
ID’ Avanzo et alj|2012j) . We note that the micro-physical param¬ 
eter values that reproduce the optical and radio fluxes underesti¬ 
mate the X-ray flux by more than one order of magnitude. This 
is not unexpected since it has already been discussed in the lit¬ 
erature that the X-ray emission of GRBs could be dominated at 
early times (typically up to half a day after the explosion) by an 
extra component which is apparently un related to the standard 
afterglow forward sh ock emission (e.g. iGhisellini et ail [2009: 
ID’Avanzo et al.ll2012l) . 

Therefore, we assume the microphysical parameters values 
that reproduce the optical and radio flux distribution of a flux 
limited sample of real bursts. 

3. Results 

Following G14 we consider the time when the OA flux reaches 
its peak (Eq. 3 in G14) and calculate the OA flux at this time 
of the peak. The light curve of the OA will start to rise when 
T — 1/ sin(0 v ; cw - 0j et ) and will peak around the time when all the 
jet is visible, i.e. T — 1 / sin(0 v iew + #jet)- The OA flux after this 
time follows the same light curve that an observer would see if 
with a line of sight within the jet opening angle. Therefore, the 
peak time corresponds to the maximum observed flux from the 
OA. According to our simulations, the time when the OA peaks 
(which depends on the burst parameters and on the viewing an¬ 
gle) has a broad distribution with a typical value of few hundred 
days after the GRB (note that this reference is purely theoretical, 
since the GRB start time is missed in the real case of an OA). 

3.1. Orphan afterglows flux density 

The cumulative peak flux density distributions of OA are shown 
in t ig- El The X-ray flux density is computed at 3 keV where the 
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Fig. 2. Cumulative flux density distribution of orphan afterglows at three characteristic frequencies: R band for the optical (red line), 3 keV for 
the X-ray (blue line) and 443 GHz for the mm band (green line - representative of ALMA frequency range). Upper limits of past searches of OA 
in the optical band are shown by the red symbols (to be compared with the red solid line). The optical R flux distribution is representative of the 
OA at z < 4.5 because, at higher redshift, their R band flux suffers from Lyo- suppression. An optical extinction (according to the distribution of 
Covino et al. 2013) has been applied to the optical fluxes. 


tending, in less than 10% of the cases, to values larger than 1-2 
mag. We note that this assumption is based on the Ay measured 
in bursts observed within their jet opening angle (i.e. close to 
the jet axis). Eventually, dust destruction by the GRB X-ray/UV 
flash could reduce the optical absorption close to the line of sight 
dPerna et al.ll2003h . OA, observed at large viewing angles, could 
therefore have a larger Ay than what we assume. Overall, this 
effect would further reduce the optical fluxes of OA observed at 
large viewing angles and consequently reduce the predicted rates 
of OA detectable in optical surveys. For the optical we have con¬ 
sidered only OA at z < 4.5 since we expect that the optical emis¬ 
sion is fully absorbed by the Ly a absorption at larger redshifts. 
Also shown in Fig. [2] are the range of variation of the flux density 
distributions obtained by varying the micro-physical parameters 
(as described in §2). 

All the flux density distributions shown in Fig. [2] represent 
the flux at the time when the OA reaches its maximum emission. 
On average, OA should reach its maximum flux hundreds of days 
after the burst. At these times the afterglow emission spectrum 
peaks at relatively low frequencies (in the mm and radio band). 
This accounts for the relative normalizations of the flux cumu¬ 
lative distributions at different frequencies in Fig. [2] The three 
curves, however, converge at very low fluxes (not shown in Fig. 
[2] for clarity) to the total rate of OA which is set by the nor¬ 
malization of the population synthesis code to the rate of GRBs 
detected in the y-ray band by Swift, BATSE and Fermi (G13). 
This can be already seen in the 443 GHz curve of Fig. [2] and 
corresponds to 10 OA deg~ 2 vr 1 . 



Log(F,J [mJy] 

Fig. 3. Average duration of the simulated population of OA with flux 
above the corresponding x-axis value. The bars represent the lcr scatter 
around the average. Typical timescales are shown by the dashed hori¬ 
zontal lines (as labelled). Linear fits are shown by the dotted lines (fit 
parameters are reported in Tab.Q] 


photoelectric absorption by metals in the Galaxy and host galaxy 
is negligible dCamnana et al.l[2012 ). For the optical R flux den¬ 
sity, we assumed an Ay according to the distribution obtained 
by the analysis of the BAT6 sample (ICovino et al.ll2013l) . This 
is an asymmetric distribution of Ay peaking < 0.5 mag and ex¬ 



F [mJy] 
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Table 2. Transient surveys in the optical and X-ray bands. On-going and future surveys are marked in boldface. Parameters of the optical surveys, 
field of view (FOV), cadence, limiting flux Fu m , coverage and lifetime are from the compilation of lRau et alJ 1 20091) . The rate of orphan afterglow 
Roa above the survey limiting flux is obtained through the flux density distribution reported in Fig. [2] The average OA duration above this flux 
limit (T) is derived from Fig.[3]and from the parameters of the linear fits reported in Tab.[j](the minimum and maximum durations are shown 
in square parentheses). The last column shows the number of O A per y ear detectable by the reported surveys. For the X-ray the sky coverage is 
intended for 24 h. * see http://www.ptf.caltech.edu/ztf and lBellin (2014 ). 


Survey 

FOV 

Cadence 

F Um 

Coverage 

Lifetime 

Roa 

(T) 

# OA 


(deg 2 ) 


(mly) 

(deg 2 night -1 ) 

days 

(deg” 2 yr 1 ) 

days 

yr -1 

PTF 

7.8 

lm-5d 

1.17xl0~ 2 

1000 


1.5xl0- 3 

1[0.2-3.8] 

1.5 

ROTSE-II 

3.4 

Id 

1.17x10 1 

450 


5.2xl0- 4 

0.4[0.1-1.7] 

0.1 

CIDA-QUEST 

5.4 

2d-lyr 

4.60xl0~ 2 

276 


8.0xl0- 4 

0.5[0.1-2.3] 

0.1 

Palomar-Quest 

9.4 

0.5h-ld 

1.17xl0" 2 

500 

2003-2008 

1.5xl0- 3 

1 [0.2-3.8] 

0.8 

SDSS-II SS 

1.5 

2d 

2.68xl0~ 3 

150 

2005-2008 

3.2xl0- 3 

1.6[0.4-6.3] 

0.8 

Catilina 

2.5 

10m-lyr 

4.60xl0~ 2 

1200 


8.0xl0- 4 

0.6[0.1-2.4] 

0.6 

SLS 

1.0 

3d-5yr 

5.60xl0~ 4 

2 

2003-2008 

5.2xl0- 3 

2.8[0.8-ll] 

0.03 

SkyMapper 

5.7 

0.2d-lyr 

7.39xl0~ 2 

1000 

2009-... 

6.4xl0- 4 

0.5[0.2-2.0] 

0.3 

Pan-STARRSl 

7.0 

3d 

7.39xl0~ 3 

6000 

2009-... 

2.0xl0- 3 

1 [0.3-4.4] 

12 

LSST 

9.6 

3d 

4.66xl0~ 4 

3300 

2022-... 

5.1xl0- 3 

3[0.8-ll] 

50 

Gaia 

0.5x2 

20d 

3.00xl0- 2 

2000 

2014-2019 

io - 3 

1[0.5-5] 

2 

ZTF * 

42.0 

Id 

2.00xl0~ 2 

22500 

2017-... 

l.lxlO- 3 

0.8[0.4-4.8] 

20 

RASS 

3.1 


4.00xl0~ 3 

12000 

6 months 

8.0x10 4 

1 [0.3-4.4] 

10 

eROSITA 

0.8 

6 months 

2.00xl0~ 6 

4320* 

4 years 

3.0xl0- 3 

2[0.5-6.5] 

26 


3.2. Orphan afterglows timescales 


Surveys can detect orphan afterglows as transient events when 
the OA flux is above the survey limiting flux. The latter, there¬ 
fore, determines the rate of detectable OA. However, the sur¬ 
vey limiting flux will also determine the OA characteristic du¬ 
ration ( T ). At a fixed frequency, the duration of OA above the 
survey limiting flux is longer the deeper is the survey. We de¬ 
fine (T) as the time interval during which the OA flux is larger 
than the survey flux limit Fu m . Fig. [3] shows the average dura¬ 
tion of OA above Fi; m for the three frequencies considered as 
reference in this work. In general, given the typical flux limits 
of optical and X-ray surveys (see also §4), OA will appear as 
daily transients. At GHz frequencies, instead they will be much 
slower transients with duration of even tens-hundreds of days 
(see also [Ghirlanda et akll2014l) . Note that this timescale repre¬ 
sents the duration of the OA above a given survey limit and it is 
only due to the instrumental limit. This should not be confused 
with the timescale of the peak of the OA with respect to the GRB 
which is due to the combination of the burst geometry (opening 
angle and viewing angle) and the hydrodynamics (i. e. the decel¬ 
eration of the fireball, i.e. mainly set by its kinetic energy, initial 
bulk Lorentz factor and circumburst density). In general, the dis¬ 
tribution of the time of the peak of the OA is centred around a 
few hundreds days. However, at these times the optical emission 
of OA is extremely faint, so that any conceivable optical survey 
(also the deepest which will be performed in the future) will de¬ 
tect those OA which peak at relatively early times, i.e. between 
1 and 10 days after the trigger, which are the brightest within the 
population. 

Fig. [3 shows that, at any frequency, the OA duration (T) in¬ 
creases as the survey limit deepens (i.e. with decreasing survey 
limiting flux Fn m ). Tab. Q] shows the parameters (slope m and 
normalization q ) of the linear fit (dotted lines in Fig. [3]) to the 
data shown in Fig.[3]for the three characteristic frequencies. 


4. Orphan afterglows detection rate 

In this section we compare our results with past searches for 
OA in the optical and X-ray band and show specific predic¬ 
tions for on-going or planned surveys in these bands. We also 
consider forthcoming large projects like the Large Synoptic Sky 
Telescope (LSST) and the extended ROentgen Survey with the 
Imaging Telescope Array (eROSITA) which will conduct almost 
all sky surveys in the optical and X-ray bands, respectively. For 
the radio band. lGhirlanda et al .l J20141) showed that the OA rates 
are consistent with the (upper) limits of past radio surveys which 
did not detect any credible orphan afterglow. Forthcoming radio 
surveys like the VAST/ASKAP at 1.4 GHz or the MeerKAT or 
EVLA at 8.4 GHz could detect 3xl(T 3 and 3x10“' OA deg- 2 
yr 1 , respectively. The deeper SKA survey, reaching the yuJy flux 
limit, could detect up to 0.2-1.5 OA deg~ 2 yr 1 ( [Ghirlanda et alJ 
|2014|) . Here we report the predictions for the optical and X-ray 
surveys. 

4.1. Optical surveys 

Among past searches for orphan afterglows in the optical, 
iRvkoff et al.l (l2005l> used the Robotic Optical Transients Search 
Experiment III (ROTSE-III). Over a period of 1.5 year, they 
identified no credible GRB afterglow. They place a 95% upper 
limit on the OA rate of 1.9 deg 2 yr -1 at R = 20. The Deep Lens- 
ing Survey (DLS - iBecker et al.l2004l) provides a (less constrain¬ 
ing) limit of 5.2 deg~ 2 day' 1 for transients with typical duration 
of a few ksec and 19.5< R <23.4. iMalacrino et al.l d2007l) ob¬ 
tained a more stringent upper limit from the CFHTLS Very Wide 
survey: excluding that the three transient they find are GRBs 
(iMalacrino et alJl2007l) . an upper limit of 0.24 deg~ 2 yr -1 down 
to R = 23 can be placed. The ROTSE-III and CFHTLS limits 
are shown in Fig. [2] (filled red symbols) and they are consistent 
with the rate for the optical band predicted by our model (solid 
red line in Fig. [2]). Also, no credible OA was_ found in the Faint 
Sky Variability Survey project (Ivreeswiikll2002h . 
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The present and future major surveys in the optical are shown 
in Tab. [2] Most of the optical survey parameters are obtained 
from [Rau et all d2009t) . In Tab. [2] we report the survey name 
(Col.l), its field of view (FOV) and its cadence (Col. 2 and 3). 
The limiting flux density and the coverage, representing the sky 
area covered per night, are reported in Col. 3 and Col. 4, respec¬ 
tively. Through Fig.0we can derive the rate Rqa of OA that have 
their peak flux density above each survey limiting flux. This is 
reported in units of deg -2 yr in Col. 7 of Tab. [2] Similarly, from 
Fig. [3] it is possible to estimate the average duration (T) of the 
OA above the survey limiting flux (Col. 8 in Tab. [2}. In square 
parentheses we indicate the upper and lower estimates of the av¬ 
erage duration (i.e. corresponding to the lcr error bars in Fig. [3}. 
We derive the rate of OA (expressed in number of OA per year - 
last column of Tab. [2]) in a given survey as Noa = Roa x C x (T), 
where C is the fraction of the sky covered by the survey per night 
(coverage in Tab. [2]) 

4.2. X-ray surveys 

Searching for GRB afterg lows in X-ray surveys led to the dis - 
covery of few flare stars (iGrindlavl 1 1 9991 : IGreiner et al. 2000). 
The 27 X-ray transients, detected in the 5.5 year survey of Ariel 
V (iPve & McHardv! 1 9831) . provide a conservative upper limit 
of 1.15 x 10_ 3 4 deg -2 yr -1 corresponding to a flux ~0.06 mJy 
dGrindlavll 1999 "). This is consistent with our predictions for the 
X-ray band (solid blue line in Fig. [2} at the same flux limit. 

Among the widest X-ray surveys, the ROSAT All-Sky Sur¬ 
vey covered the full sky reaching a limiting flux of 4x10 M erg 
cm' 2 s -1 (0.5-2 keV) in almost 6 months. This flux limit (assum¬ 
ing a spectrum with photon index -2) corresponds to a flux den¬ 
sity of ~ 4 x 10 -5 mJy which, according to our estimates, gives a 
rate Roa ~ 8 x 10 4 deg -2 yr' 1 . The RASS scan procedure cov¬ 
ered a full sky circle of width 2 deg every orbit corresponding 
to ~ 12000 deg 2 day -1 . According to our estimates (Fig. [2] and 
Tab.0, the typical duration of the OA above the RASS flux limit 
should be of ~ 1 day so that the expected OA number should be 
~ 4.8 during the 6 months survey lifetime. This result is consis¬ 
tent with the estimate of iGreiner et al.l (j200C)t ). They effectively 
searched in the RASS for GRB afterglows and concluded that of 
the 23 candidates only a few could be due to GRBs, most of the 
others being flare stars. The second release of the RASS, 2RX50 
has been extended to a flux limit a factor of 4 deeper than the 
first release. Therefore, we expect to have ~12 OA in the 2RXS. 

In the X-ray band Chandra and XMM-Newton have per¬ 
formed deep surveys but, due to their small field of view, at 
the expense of a relatively s mall portion of the sky explored 
(see [Brandt & Hasingerll2005l) . The observing strategy in these 
surveys was not a scanning mode as the RASS, but rather the 
combination of pointed repeated observations of the same field. 
Therefore, it is difficult to reconstruct the overall sky coverage. 
As a gross estimate we can compute the expected number of 
OA by multiplying the predicted rate (according to our results 
of Fig. 0 times the area of the sky covered. We stress that this 
is an overestimate of the number of OA that could be detected 
by these surveys. Among the deepest surveys, the 2Ms Chan¬ 
dra Deep Field North covered 0.13 deg 2 in the 0.5-8.0 keV band 
down to a flux limit of ~ 10 16 erg cm -2 s 1 (Alexander et al J 
l2003i ). To such a flux limit we predict less than lO -2 OA yr -1 . 
Similar rates are expec ted in the XMM-Newton Large Scale Sur¬ 
vey dPierre et al.ll2004l ) which, with a sensitivity of ~ 5 x 10 -15 


2 http://xmm.esac.esa.int/external/xmm_science/workshops/ 
2014symposium 


erg cm -2 s 1 (0.5-2 keV) and a 10 deg 2 of sky coverage, should 
detect at most 0.1 OA yr -1 . Both Chandra and XMM-Newton 
have performed several other surveys, but, despite with larger 
sky coverage than those mentioned above, at the expense of their 
sensitivity (e.g. the XMM-Newton Bright Serendipitous Survey 
- [Della Ceca et al.l (120041) 1. According to our predictions, there¬ 
fore, there seems to be no chance to have detected any orphan 
afterglow in current deep X-ray survey^ 

4.3. Orphan Afterglow distinguishing properties 

The most important question to address is how to distinguish 
OA from other transients when they will be detected in large 
sky surveys. As shown in Fig. 0OA will appear as daily tran¬ 
sients in optical and X-ray surveys, given the typical flux limits 
of current and forthcoming surveys (see Tab. 0, and many other 
extragalactic sources will have similar duration. They will show 
a decaying light curve with a temporal slope which is that of typ¬ 
ical GRB afterglows, i.e. oc r s with 8 -1-2 but not uniquely 
characterising GRBs as a class. The lack of any associated high 
energy y-ray counterpart will hamper their classification as or¬ 
phan afterglows of GRBs. The success of a transient survey is 
that of classifying, after discovery, the detected transients and to 
this aim a dedicated follow up program is fundamental. 

A first way to assess the orphan afterglow nature of the tran¬ 
sient is by a systematic optical photometric and spectroscopic 
follow-up. The optical/X-ray light-curves, especially if the OA 
that is detected is still quite bright and/or before its light-curve 
peak, can be a very useful tool for a preliminary source classifi¬ 
cation. Indeed, its shape and its decay power-law index should be 
different from those of SNe or blazars. In addition to the ground- 
based optical facilities, the future satellite SVOM dBasa et al.l 
120081) . to be operational at the same epoch of LSST, will be able 
to perform simultaneous X-ray and optical observations of OA 
candidates. 

A final identification is given by optical spectroscopy. The 
spectral continuum and the absorption lines present in the 
optical-NIR spectra of the afterglow are very different from 
those typically found in SNe or blazars dFvnbo et al.l 2009; 
IChristensen et al J1201lh . Nonetheless, spectra with a sufficient 
signal-to-noise ratio are needed. To date this requirement is ful¬ 
filled down to R ~22 with a reasonable amount of integration 
time (~2hrs) (e.g.: with X-shooter at ESO/VLT). Future larger 
telescope will make possible to obtain similar results also for 
fainter objects. Therefore, with the possible increase of transient 
detections from future surveys, a considerable amount of tele¬ 
scope time may be necessary for a systematic optical follow-up 
of unknown transients and of potential OA candidates. Among 
projects which could substantially contribute to the broad band 
(optical to NIR) spectroscopy of transients discovered in surveys 
there is the Son Of X-Shooter (SOXS - RI. Campana) which has 
been proposed for the NTT with a considerable number of dedi¬ 
cated nights per year and the Nordic Optical Telescope Transient 
Explorer (NTEjj. Another way to discriminate between OA and 
other transients comes from the analysis of the broad Spectral 
Energy Distribution (SED). Here we compare the typical SED 
of OA with that of potential competitors extragalactic sources 
like SNe and blazars. 


3 In the XMM/EPIC database, wh ose vari ability richness will be fully 
explored by the EXTraS project i lDe Luca et all20f5il . we do not expect 
orphan afterglows to be present. 

4 http://dark.nbi.ku.dk/news/2014/nte_is_a_go/ 
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Frequency [Hz] 

Fig. 4. Spectral energy distribution of the OA that can be detected by the LSST (pink filled region). The SED of the low power blazar BLLac (open 
circles), of the FSRQ 3C454.3 (asterisks) and of two supernovae SN1978K (open squares) and the GRB-SN associated GRB980425/SN1998bw 
(open stars) are shown by different symbols. The solid lines provide an interpolation of the data points and do not represent any physical model. 
For the two blazars we also show (dashed grey line for 3C454.3 and dashed orange line for BLLac) how their spectra would appear if they 
were at z = 2 (typical of GRBs). The (5<x) limits for a 12 h continuum observation with the SKA is shown by the yellow shaded region. 
The green shaded region marks the limiting flux of an ALMA observation (32 antennas of 12 m for 3 h of observation in dual mode - from 
https://almascience.eso.org/proposing/sensitivity-calculator). The LSST liming flux (see TabO is shown by the red square symbol. 


Since the most promising detections will be with the forth¬ 
coming LSST (Ivezic et al. 2008 - see Tab.0, we have consid¬ 
ered only the OA that will be detected by this survey. We predict 
a rate ~ 50 OA per year. The overall SED (i.e. the convolution 
of the SEDs of all OA detectable by the LSST survey) is shown 
by the hatched pink region in Fig. [4] The typical SED of OA de¬ 
tectable by LSST peaks in the 10 11-13 Hz range. The spectrum 
below the peak, in the GHz down to the MHz range scales cc y 2 . 

Possible extragalactic variable sources that could compete 
with GRB orphan afterglows in brightness, frequency of discov¬ 
ery and timescales are supernovae and blazars. Fig. [4] shows the 
SED of two blazars: the Flat Spectrum Radio Quasar (FSRQ) 
3C454.3 and BLLac itself as representative of the respective 
classes (see iGhisellini et al.l l2010i) . Two supernovae are also 
shown: SN 1978K as a possible representative of highly lumi¬ 
nous supernovae and the GRB980425/SN1998bw (iGalama et alJ 
Il998h for the class of GRB-SNe associated. For all these 
sources, we report their SED as obtained by multifrequency ob¬ 
servational campaigns and retrieved from Italian Space Agency 
(ASI) Science Data Center Sed Builder too0. The solid curves in 
Fig- E] are not physical models but only illustrative of the overall 
broad band spectral energy distribution of these classes of ob¬ 
jects. For the blazars we also show how their SED would be like 
if they were shifted at z = 2, i.e. at the typical distance of long 
GRBs. 


5 http://tools.asdc.asi.it 


For comparison in Fig. [4] we show the LSST flux limit (red 
square symbol). The OA that can be detected by LSST when 
their jet emission is fully visible by the off-axis observer will 
have their peak frequency already below the optical band, in 
the mm region. This is because the peak of the OA emission 
is reached several months after the burst (§3 - see also G14). 
Furthermore, Fig. 0] shows that their emission in the MHz/GHz 
region is still in the self absorbed regime. Differently, BLLacs 
and SN emission like 1998bw or 1978K are characterised by a 
softer spectrum in the radio band than the typical OA detected 
in an optical survey like the LSST. Therefore, the follow up of 
these transients in the mm and GHz bands will characterise their 
different SED. 

5. Discussion 

Among previous wo rks in the literature which e stimated the de¬ 
tection rate of OA, iTotani & Panaitescui (12002 ) considered 10 
GRB of the pre-Swift era with well monitored afterglow light 
curves as templates. By assuming different off-axis viewing an¬ 
gles they estimated the rate of OA in the X-ray, optical and radio 
band. Their predictions were based on a very small number of af¬ 
terglows mostly representative of the bright afterglow population 
of GRBs. Similarlv lZou et all d2007l) adopted a set of fixed phys¬ 
ical parameters (kinetic energy and micro-physical parameters) 
and allowed only for a possible distribution of jt. They predict a 
rate of 1.3x10 2 deg -2 yr“ 1 for OA brighter than R - 20 which 
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is consistent with our findings at the corresponding flux (red 
solid curve in Fig. [2]). However, their flux distributions appear 
steeper than our model at low fluxes, thus predicting a larger rate 
of OA in surveys that go deeper than the above limit. We con¬ 
sider that our estimates better represent the low flux end of the 
OA distribution since our code is calibrated with the entire GRB 
prompt emission flux distribution and includes a more represen¬ 
tative sample of afterglows to fix the micro-physical parameters. 

What is new in our model is that we predict the properties of 
OA based on the observed properties of GRBs in the y -ray band 
considering as constraints the flux and fluence distribution of the 
population of GRBs detected by Swift, BATSE and Fermi. Since 
the y-ray energy detected in the prompt emission is a proxy of 
the kinetic energy driving the afterglow deceleration, our sim¬ 
ulated population of bursts includes both high and low kinetic 
power bursts. The choice to reproduce the afterglow flux dis¬ 
tribution of the complete BAT6 sample of Swift bursts, despite 
being composed by relatively bright events, ensures that we are 
extending the flux distribution of the synthetic GRB population 
to the low end better than what could be done with the limited 
number of GRBs detected in the pre-Swift era. 

We have assumed that GRBs have a jet with a top-hat uni¬ 
form structure, i.e. the kinetic energy and the bulk Lorentz 
factor are constant within the jet opening a ngle. Alternatively 
dRossi et al.ll2002llZhang & Meszarosll2002l) GRB jets could be 
structured, i.e. the kinetic energy (and possibly also) the bulk 
Lorentz factor depend (as a power-law or exponentially) on the 
angle from the jet axis. In the former scenario, considered in 
this work, orphan afterglows are naturally expected to dominate 
the number of GRBs in the Universe (considering a typical jet 
opening angle of few degrees). Even if the top-hat uniform jet 
emission can be seen when 0 v iew^f>. its flux decreases drastically 
for off-axis observers so that we can consider valid the approx¬ 
imation so far understated that we do see the prompt emission 
only of GRBs whose uniform jet is pointed towards the Earth 
(i.e. when (ficw<b)- I' 1 the structured jet model, instead, there is 
always a portion of the jet that is pointing to the observer. There¬ 
fore, the observed GRB properties depend only on the viewing 
angle 0 v j ew so that orphan afterglows should not exist in princi¬ 
ple, since even at large angles from the jet there is jet emission 
that can be seen dSalafia et al.ll2015h . However, also in this sce¬ 
nario OA could still be present if either the prompt emission at 
large angles is below any detector threshold or if the jet is uni¬ 
form within a relatively narrow core and highly structured (i.e. 
with a steeply decreasing energy profile) outside it, as suggested 
by recent results from the modelling of the luminosity function 
dPescalli et al.l [20151) . In this model the detection rates of OA 
would be, anyw ay, sm aller than in the top-hat model adopted 
here as shown in lRossi et al. ( 2008 ). 

Finally, note that in Fig. Q] we have shown that our after¬ 
glow model does not reproduce the X-ray early flux observed 
in the BAT6 sample of GRBs. This is due to the fact that at 
early times the X-ray emission has been shown to be incon¬ 
sistent with (i.e. larger than) the forward external shock model 
(e.g. IGhisellini et alJl2009l) . However, we have used the after¬ 
glow emission of the population of OA (as shown in Fig. [2} in 
predicting the rate of detection of OA by forthcoming X-ray sur¬ 
veys. Indeed, the time when the brightest OA peak in the X-ray 
(i.e. those that will likely be above the threshold of forthcom¬ 
ing X-ray surveys) is a few days after the prompt emission i.e. 
when also in the X-ray th e emission turn out to be dominated by 
the afterglow component (IGhisellini et al.l2009tlD’ Avanzo et alj 
1201 21 ) . 


6. Conclusions 


We computed the emission properties of the population of or¬ 
phan afterglows in the optical and X-ray band. Our simula¬ 
tion procedure relies on what we have so far observed for a 
well defined and complete sample of GRBs detected by Swift 


emission have been extensively studied (Campana et al. 

2012 

D’Avanzo et al] 

2012S: Covino et al. 20131; Melandri et al 

2014 

Ghirlanda et al. 

2013a). 


Our model allows us to predict the expected rate of detection 
of OA in past, current and future optical and X-ray_surveys_(Tab. 
13, Fo r a simila r work for the radio band see iGhirlanda et al.l 
120141) (see also lMetzeer et al.l (120151) ). Most past and on-going 
optical surveys have small chances to detect OA. Among these, 
the Palomar Transient Factory (PTF - lLaw et akl (12009!) ) could 
marginally see one OA per year (consistent with Rau et all2 009. 
predictions) given its relatively low sensitivity compensated by 
the large portion of the sky covered per night (10 3 deg 2 ). In¬ 
stead, according to our model, an optical survey like that of 
Pan-STARRS 1 which will cover 6000 deg 2 per night could al¬ 
ready detect a dozen of OA per year. Larger detection rates are 
expected for the forthcoming development of the PTF survey. 
The Zwicky Transient Facility (lBellmll2014l) . which is designed 
specifically for transients discovery, will cover about 22500 
deg~ 2 per night down to a limiting magnitude ~20.5. We expect 
that it will detect ~20 OA yr 1 . A considerably larger number 
of OA will be accessible by the Large Synoptic Sky Telescope 
survey (LSST - llvezic et al.l200 8). The telescope will have a 9.6 
deg 2 field of view and will be able to survey 10 4 deg 2 of the sky 
every three nights down to a limiting magnitude for point sources 
R ~24.5. With these parameters we estimate it could detect 50 
OA yr _1 . 

An interesting prediction concerns the Gaia satellite 
(tLindegrenll2010l) . It will carry two telescopes each one with a 
field of view of 0.7°x0.7° and will scan an angle of 360° ev¬ 
ery six hours. Therefore, it will cover ~2000 deg 2 per day per¬ 
forming a survey down to a limiting flux of 0.03 mJy. Accord¬ 
ing to our model Roa ~ 10~ 3 deg~ 2 yr~* at this flux limit so 
that we predict that Gaia will detect about 10-15 OA in its 5 
year mission. This estim ate is consistent with that reported in 
llaneli & Gombod(l2011l) . 

Given the depth of forthcoming optical surveys we expect 
that OA will have a typical redshift z ~2. At such distances the 
typical GRB host gala xy s hould be fainter than the LSST limit¬ 
ing magnitude dHiorth et al.112012l) . 

The difficulty will be to disentangle these OA from other 
(galactic and extragalactic) transients that will be detected with 
similar flux and temporal behaviour. The follow up in the opti¬ 
cal and X-ray will secure the sampling of the light curve which 
could be the first hint to the OA nature (with respect to potential 
other extragalactic transients like supernovae and blazars). The 
availability of dedicated facilities or assigned observing time at 
different ground based facilities will be crucial in this respect. 
Optical/NIR spectroscopy will discriminate extragalactic tran¬ 
sients (e.g. lGorosabel et al.ll2002l) . low frequency (mm and radio 
bands) observation^ could be used (as shown in Fig. S to dis¬ 
tinguish among possible competing transients sources. For par¬ 
ticularly low redshift transients, the search for the host galaxy 
could also provide further clues on their nature. 

Here we have computed the OA flux cumulative distribution 
at the reference frequency of 443 GHz which is one of the ffe- 


6 The LSST will start operating approximately in the same period of 
the Square Kilometer Array (SKA). 
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quencies covered by e.g. ALMA. We have verified that the few 
hundred GHz range is where OA are brightest considering the 
typical timescales when they become visible (§2). Indeed, at 40 
GHz and 4000 GHz the flux cumulative curves in Fig. [2] lie be¬ 
low the solid green curve representing the flux density at 443 
GHz. The Herschel/SPIRE survey ATLAS0, one of the widest 
covering a total of 500 deg 2 , is limited by the confusion limit 
of 5-7 mJy at 250-500 pm so that we expect less than 0.1 OA 
yr _I . Spitzer SW1RL0 observed six fields in the northern and 
southern sky with typical areas between ~4.2 deg 2 and ~ 12 deg 2 
with higher sensitivities of few tens of pjy in the low frequencies 
channels (IRAC) at 3.6pm and 4.5pm. These fields were covered 
on different timescales between 1 and 6 days. According to our 
model we expect a rate of less than one OA per year in such 
fields above the deepest flux limits of this survey. 

Among forthcoming X-ray surveys we consider the ex¬ 
tended ROentgen Survey with the Imaging Telescope Array 
(eROSITA - Merloni et al. 2012) which will cover the full sky 
up to 10 keV with a flux limit ~ 2 x 10~ 14 erg -1 cm” 2 s” 1 in the 
0.5-2 keV band. Therefore, ~ 3 x 10 ! deg” 2 yr” 1 OA should 
be reachable by this survey (Tab. [2}. According to the planned 
scanning strategy, a full circle of width 2 degree will be covered 
every four hours. This corresponds to ~ 4320 deg 2 d ay” 1 . The 
expec ted OA number is ~26 yr” 1 (but see also lKhabibullin et all 
120121) . A larger number o f OA (by a factor 2 ) could be reached 
by the WFXT survey (e.g. lRosati et alJl201 ll) . 
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